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ABSTRACT 

Context. Models of stellar structure and evolution can be constrained by measuring accurate parameters of detached eclipsing binaries 
in open clusters. Multiple binary stars provide the means to determine helium abundances in these old stellar systems, and in turn, to 
improve age estimates. 

Aims. Earlier measurements of the masses and radii of the detached eclipsing binary V20 in the open cluster NGC 6791 were ac- 
curate enough to demonstrate that there are significant differences between current stellar models. Here we improve on those results 
and add measurements of two additional detached eclipsing binaries, the cluster members V18 and V80. The enlarged sample sets 
much tighter constraints on the properties of stellar models than has hitherto been possible, thereby improving both the accuracy and 
precision of the cluster age. 

Methods. We employed (i) high-resolution UVES spectroscopy of VI 8, V20 and V80 to determine their spectroscopic effective tem- 
peratures, [Fe/H] values, and spectroscopic orbital elements, and (ii) time-series photometry from the Nordic Optical Telescope to 
obtain the photometric elements. 

Results. The masses and radii of the V18 and V20 components are found to high accuracy, with errors on the masses in the range 
0.27-0.36% and errors on the radii in the range 0.61-0.92%. V80 is found to be magnetically active, and more observations are needed 
to determine its parameters accurately. The metallicity of NGC 6791 is measured from disentangled spectra of the binaries and a few 
single stars to be [Fe/H]= +0.29 ± 0.03 (random) ± 0.07 (systematic). The cluster reddening and apparent distance modulus are found 
to be E(B - V) = 0.160 ± 0.025 and (m - M)v = 13.51 ± 0.06 . A first model comparison shows that we can constrain the helium 
content of the NGC 6791 stars, and thus reach a more accurate age than previously possible. It may be possible to constrain additional 
parameters, in particular the C, N, and O abundances. This will be investigated in paper II. 

Conclusions. Using multiple, detached eclipsing binaries for determining stellar cluster ages, it is now possible to constrain param- 
eters of stellar models, notably the helium content, which were previously out of reach. By observing a suitable number of detached 
eclipsing binaries in several open clusters, it will be possible to calibrate the age-scale and the helium enrichment parameter A Y/AZ, 
and provide firm constraints that stellar models must reproduce. 

Key words. Open clusters: individual NGC 6791 - Stars: evolution - Stars: binaries: spectroscopic - Stars: binaries: eclipsing - 
Techniques: spectroscopy - Techniques: photometry 



> 1. Introduction 



The open cluster NGC 6791 is interesting for several reasons. 
The main scientific motivation for studying this cluster is that it 
is one of the oldest an d, at the same time , the most metal-rich 
open clusters know n (Origli aet al.1 l2006t ICarretta et al.l [2007: 
Anth onv-Twarog et al.l 120071) . and thus very important for un- 
derstanding chemical evolution. In addition, it is a particularly 
populated open cluster, with stars in all st ages of evolution 
from the main sequence to the white dwarfs dKing et al.l l2005: 



Bedin et al. 2005L 2008 


), as well as with numerous variable 


stars (Bruntt et al. 2003; 


Mocheiska et al.ll2002r.lde Marchi et al.l 


2007). 



Send offprint requests to: K. Brogaard, e-mail: kfb@phys.au.dk 

* Based on observations carried out at the Nordic Optical Telescope 
at La Palma and ESO's VLT/UVES ESO, Paranal, Chile (75.D-0206A, 
77.D-0827A, 081.D-0091). 

** Tables 11-22 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



Despite a number of studies that have presented well- 
calibrated photo metry and high-precision colour -magnitude di- 
agrams (CMDs) dKing et al.ll2005l IStetson et al.ll2003l hereafter 
SBG03), the age of NGC 6791 has remained very uncertain until 
recently because of correlated uncertainties in distance, redden- 
ing, colour-temperature transformations, and metallicity. 

It is widely appreciated that detached eclipsing binaries offer 
the possibility of determining accurate (and precise) masses and 
radii for the system components, nearly in dependent of model 
assumptions dAndersenlll991l : iTorres et all 1201 Oh . If the binary 
resides in a star cluster, and one or both of its components are 
close to the turn-off, it is possible to put tight constraints on 
the age of the system by comparing the position of the primary 
and secondary in a mass-radius (MR) diagram to theoretical 
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Table 1. Names and coordinates for the eclipsing binaries. 



Name Alternative name RA(2000.0) DEC(2000.0) 



V18 V565 Lyr 19 : 20 : 49.38 +37 : 46 : 09.3 

V20 V568Lyr 19:20:54.30 +37:45:34.7 17.5 

V80 - 19 : 21 : 06.48 +37 : 47 : 27.8 



isochrones. For stellar clusters, such an analysis has some sig- 
nificant advantages: the determination of the masses and radii 
is independent of the usual uncertainties such as reddening and 
distance. Furthermore, since the comparison to models is car- 
ried out in the MR diagram, one avoids the difficult process of 
transforming the effective temperatures and luminosities of the 
models to observed colours and magnitudes. Thus, determining 
cluster ages in the MR diagram allows a direct confrontation be- 
tw een observations and theory. 

iGrundahl etal] (2008; hereafter GCH08) showed that us- 
ing this method with their measurements of the cluster mem- 
ber eclipsing binary V20, they could determine a precise cluster 
age with an error of only +0.3 Gyr for a given stellar model. 
However, they were unable to determine which of the models (if 
any) to trust, because the difference in predicted age due to the 
specific stellar model adopted was about four times greater than 
their measurement precision. 

Here we undertake an analysis of three detached eclipsing 
binary systems, VI 8, V20 and V80, in NGC 6791, and deter- 
mine accurate masses and radii for the components of two of 
the systems. We also measure spectroscopic T e ff and [Fe/H] val- 
ues from disentangled spectra of the binary stars. This is used, 
together with the cluster CMD, to demonstrate how multiple 
cluster member eclipsing binary measurements constrain stel- 
lar models and cluster parameters, like age and helium content, 
better than previously possible. 

2. Targets 

The detached eclipsing binaries analysed here are VI 8, V20 
and V80, with names, alternative identifications, and coordi- 
nates given in Table [TJ along with finding charts in appendix 
A. These systems we re first identified as e clipsing binaries by 
iRucinski et all d 19961) (V18 and V20) and iBruntt et all d2003b 
(V80). Their observed locations in the (B - V), V CMD are 
shown in Fig. Q] together with calculated locations of the indi- 
vidual components as determined in Sect. 13.31 

3. Photometry 

The photometric data for V18 and V80 consists of V (Johnson) 
and R (Cousins) CCD observations from the 2.56m Nordic 
Optical Telescope (NOT) and its ALFOSC instrument. In order 
to reduce the readout time, while still including both V18 and 
V80 in each frame, we limited the readout section of the CCD to 
half its full size, thus covering a field-of-view of 6'28" x 3' 12". 
We refer to the telescope homepage for further information about 
the ALFOSC instrumenfl 

Observations were carried out on 9 nights between May and 
August 2009. Most observations were obtained on nights when 
both systems showed eclipses, in order to make efficient use of 
telescope time, resulting in a total of 580 and 538 exposures in V 
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Fig. 1. The eclipsing binaries VI 8, V20 and V80 in the CMD of 
NCG 6791. Squares are the total system light, for V20 includ- 
ing a third light {I3). Circles represent the calculated positions 
of the individual components as determined in this paper from 
light curve analysis and the assumption that all components are 
positioned on the cluster main sequence, p, s, and Z3 indicate pri- 
mary, secondary, and third light components, respectively. 

and R, respectively. For all observations we employed an expo- 
sure time of 240s in V and 180s in R. The photometric data for 
V20 used in this paper are the same as those by GCH08. That 
data set also contained an egress of the primary eclipse of VI 8, 
which was missed during the 2009 campaign, thereby allowing 
us to obtain complete coverage of both eclipses and better out- 
of-eclipse coverage in the V light curve for this object. 

The bias frames and flat fields, used in the data re- 
duction, were obtained during evening twilight on each 
observing night. All photometry was carried out with 
DAOPHOT/ALLSTAR/ALLFRAME dStetsonlll987l 11994 and 
transformed to a common coordinate system using MATCH 
and MASTER (R Stetson, private comm.). Each frame was 
processed using a point-spread function (PSF) calculated from 
about 200 of the brightest stars in the field. 

3.1. Light curves and standard indices 

With the instrumental photometry in hand, we proceeded to 
transform the observations to the V and R standard system us- 
ing the same methodology as in GCH08. SBG03 put a large ef- 
fort into the transformation of BVI magnitudes onto the stan- 
dard system, and we have used the available V photometry from 
this source as internal standard stars. For each frame, a linear 
transformation from instrumental magnitudes to standard mag- 
nitudes was calculated, using (B - I) (available from SBG03 for 
all stars in the field) as the colour term. Subsequently we av- 
eraged the coefficient for the colour term for all frames (for a 
given filter) and used this for the final determination of the ze- 
ropoint for each frame. The BVI photometry listed in Table [2] 



,V20_total 



'£ini*i' • B Vf8_tota 

aYSO-totai 

#V80p".- - ' 
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is therefore on the same system as SBG03. The accuracy of 
the photometric zeropoint is in the range (P01 to 0^02 as men- 
tioned in SBG 03. Our source for /? stan dard magnitudes was the 
photometry by Mocheiska et al. (2005), as downloaded from the 
PISCES homepageQ. Apart from the different source of standard 
magnitudes, the proced ure followed was identical to that for V. 
Mocheisk a et al.l J2005) reports an offset in their V photometry 
of 0.047 relative to SBG03, and therefore the (V - R) colours in 
Table|2]may not be very accurate. 

The V and R light curves of VI 8 and V80 are listed in 
Tables 11-14 on CDS. A few obvious outliers and systemati- 
cally deviating measurements at dusk and dawn have been re- 
moved. V light curves of V18, V20 and V80 are shown in Fig.[2|- 
[4]with phases calculated from the ephemerides given in Eqs.QJ[3] 
(Sect. E3). For V20 we reuse the light curves from GCH08. In 
the top panels, all systems are shown on the same scale for easy 
comparison. 

Throughout the paper, the component eclipsed at the deeper 
eclipse at phase 0.0 is referred to as the primary (p), and the 
other as the secondary (s) component. 

V18 and V20 are seen to be well detached with a practically 
constant light level outside eclipses, while V80 shows signs of 
magnetic activity, most likely due to its shorter orbital period. 
For V20 and V80 the secondary eclipse occurs at phase 0.50, and 
the eclipses are of equal duration, suggesting that their orbits are 
circular. V18 has a slightly eccentric orbit, as seen by the offset 
from phase 0.5 of the secondary eclipse. 



3.2. A new measurement of the third light component of V20 

For V20, a very close companion is included in the light curves, 
meaning that a significant amount of third light is present. 
GCH08 estimated the amount of third light by finding a solu- 
tion where the total light of the system is shared among three 
stars constrained to lie on the cluster main sequence. However, 
they also mention that high-resolution imaging would allow a 
direct determination of the contribution to the total light from 
the third star. We therefore identified V20 in ACS images from 
the Hubble Space Telescope (HST), measured the magnitudes 
of the light emitted by the third star, which is well separated 
from the binary, and transformed the resultant magnitudes in the 
HST filters to V and /. To accomplish the latter, we selected stars 
bearing the closest similarity to the third light that we could find; 
specifically, stars within 0.003 mag in the (ra606 - m814) colour 
index and within 0.05 mag in »?606. The four twins to the third 
light that were identified this way, in our HST photometry, were 
located in the ground-based observations of SBG03, from which 
their B, V, and / magnitudes were determined. 

The mean and rms spread of these four stars is our mea- 
surement of the third light magnitudes in these bands. We find 
B = 20.23 + 0.01, V = 19.19 + 0.01 and / = 18.12 + 0.01. Since 
we only have four stars these error estimates could be underes- 
timated and we therefore conservatively adopt an error of +0.02 
for all bands. From these and the magnitudes of the combined 
light of V20 we calculated the third light contribution to V20, 



which is 0.164 + 0.004 in B, 0.182 + 0.004 in V 



3 L p +L s +L } 

and 0.197 + 0.004 in /. These values are significantly higher and 
much more accurate than 0.146 + 0.022 in V and 0.168 + 0.025 
in / as found indirectly by GCH08. 
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Fig. 2. Phased V light curve for VI 8. The secondary eclipse is 
shifted away from phase 0.5, showing that the orbit is eccentric. 
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-0.03-0.02-0.01 0.00 0.01 0.02 0.03 0.47 0.48 0.49 0.50 0.51 0.52 
phase phase 

Fig. 3. Phased V light curve for V20 from GCH08. The sec- 
ondary eclipse is at phase 0.5 implying that the orbit is circular. 
The secondary eclipse is total. 



3.3. Standard photometry of the binary components 

Standard, out-of-eclipse, B, V, R, / photometry for the binaries is 
listed in Table [2] together with calculated individual photometry 
for the stars from their light curve solutions and, in addition for 
B - V, the assumption that all stars of each system are located 
on the cluster CMD fiducial (as done for V20 in GCH08). For 
V80 the colours are calculated using only this assumption, which 
results in L s /L p = 0.15 + 0.05 in V, since the light curve does 
not constrain the light ratio of this system. 



3.4. Periods and ephemerides 

For V20 we adopted the ephemeris of GCH08. For V18 and V80 
we produced phased diagrams adding small shifts to the spectro- 
scopic periods in order to determine the best initial estimates for 
the periods and ephemerides. We used these initial estimates to 
determine both the epoch and the period from analysis of our V 
and R light curves (see Sect. 13.5b . For each system the periods 
obtained from the two light curves agree well and are very close 
to the initial estimates. 
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Fig. 4. Phased V light curve for V80. The secondary eclipse is 
at phase 0.5, implying that the orbit is circular. Variations in the 
overall magnitude level indicate magnetic activity. 

Table 2. Standard photometry for the eclipsing binaries and 
their components. 



Name 


V 


B- V 


V-R 


V-l 


V 1 ^fofal 


17.713 


0.940 


0.505 


0.990 


V 1 "primary 


18.247 


0.916 


0.478 


(...) 


^ 1 ^secondary 


18.740 


0.978 


0.546 


(...) 


V20 Mal 


17.340 


0.930 


(...) 


0.985 


V^20 primary 


17.695 


0.876 


(...) 


0.923 


secondary 


19.874 


1.162 


(...) 


1.246 


V20,, 


19.190 


1.040 


(...) 


1.070 


VS0 tota l 


17.886 


0.932 


(...) 


1.015 


V80 primary 


18.037 


0.897 


(...) 


(...) 


V^80 [secondary 


20.100 


1.200 


(...) 


(...) 



Notes. The reason some colour indices are missing is that the we do 
not have light curves and light ratios for all systems in all colors. For 
(B-V) we have calculated the colour indices by finding solutions where 
all components are constrained to the main sequence. This is done to be 
able to put all systems on the same CMD (Fig. [TJ. 



We adopt the following linear ephemerides for all analysis in 
this paper: 



MinI(V18)= 2454651.4506 + l&imb'SZ xE 



±1 



:10 



MinI(V20)= 2453151.6061 + H^QIS xE 

±9 ±25 

Min I (V80) = 2454652.3045 + 4<?88594 x E 



-A3 



+ 16 



(1) 



(2) 



(3) 



3.5. Photometric elements 

V18 and V20 are well detached systems with no signs of out-of- 
eclipse variability. Furthermore, the components of each system 
are small relative to their separation. Thus, they are well suited 
for accurate measurements. The light curves of V80 show signs 
of magnetic activity (Fig. |4j), and more observations are needed 
before we can hope to obtain precise and accurate photometric 
elements for this system. However, since we can already get a 
good measurement of log g, which is useful for the further spec- 
troscopic study, we include preliminary dimensions of this sys- 
tem as well. 

We adopted the simple Nel son-Davis-Etzel model 
dNelson & Davis! Il972t lEtzell 1 198 lb for the light curve 
analysis. It represents the deformed stars as biaxial ellipsoids 
and applies a simple bolometric reflection model. We have 
used the corresponding JKTEBOPlI code, which is a revise d 
and extended version of the original EBOP code (lEtzell ll98 1 ). 
The Levenberg-M arquardt minimization algorithm (MRQMIN: 
iPress et al] 1 19921) is used for the least-squares optimization 
of the parameters, and the code has been extended to include 
non-linear limb darkening and adjustment of epoch and orbital 
period. We made use of important new features in the latest ver- 
sion, namely the possibility to fit the light cur ve simultaneously 
with an externally deter mined light ratio dSouthworth et alj 
12007b . a third light ratio JSouth worth! l20loh , and esinw and 
ecosto dSouthwor th et al.l l2009t) . including their errors. In 
some o f its modes, JKTEBOP performs Monte Carlo simu- 
lations ( Southworth et al.l l2004albl) and residual shift analysis 



dSouthwor th 2008), which we used to assign realistic errors to 
the photometric elements. 

The magnitude at quadrature was always included as a free 
parameter, and the phase of primary eclipse was allowed to shift 
from 0.0. In initial JKTEBOP analysis, the epoch and orbital 
period were included as free parameters and then fixed at the 
values given in Eqs.QJO A circular orbit was assumed for V20 
and V80 as found from the spectroscopic orbits. For V20, a light 
curve analysis with the eccentricity as a free parameter also re- 
sulted in a circular orbit supporting the assumption, while the 
V80 light curve is not of a high enough quality to make such 
a test. The mass ratios, q, between the components were kept 
at the spectroscopic values (Table |6j. Gravity darkening coef- 
ficients corresponding to convective atmospheres were applied, 
and the simple bolometric reflection model built into JKTEBOP 
was used. Tests showed that these assumptions have negligible 
effect on the derived photometric elements. 

The procedure we followed in order to minimize the errors 
related to limb darkening deserves some attention. It is well 
known that the linear limb-darkening law is a poor fit to both the 
observed limb darkening of the Sun and that predicted by the- 
oretical model atmospheres. We therefore adopted a non-linear 
limb-darkening law, more precisely the square root law, since 
according to van Hamme (1993) this gives the best fit to stel- 
lar atmospheres in the temperature range of these binary stars. 
Theoretical coefficients to be used with the square root law for 
a given star have been c alculated by several authors us ing ID 
stellar atmospheres (e.g. van Ha mmd[l993l : IClaretll2000l) . They 
all depend on the T e g, logg, and [Fe/H]. Furthermore, they de- 
pend on the stellar atmosphere model used for the calculation. 
To make things more complicated it has been shown that ID 
stella r atmospheres d o not reproduce the limb darkening of the 
Sun (ISing et al.ll2.008h . and therefore theoretical limb-darkening 
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coefficients will be inaccurate in any case. In order to minimize 
the errors related to these facts, the optimal approach would be 
to fit the coefficients as part of the light curve solutions, instead 
of fixing them at theoretical values. But most light curves, in- 
cluding ours, are not of sufficient quality to fit for both a lin- 
ear and non-linear c oefficient for both stars. Instead we exploit a 
fact demonstrated by South worth et al. (2007) that for the square 
root law, the linear and non-linear coefficients are highly cor- 
related. Therefore, by fixing one coefficient and fitting for the 
other, any error in the fixed coefficient will be compensated for 
by a shift in the value of the fitted coefficient. This is the ap- 
proach we adopted, since it allowed us to use a non-linear law 
while minimizing dependence on the model and the exact stellar 
parameters. 

In the text and tables with photometric solutions we use the 
following symbols: i orbital inclination; r p = R p /a; relative ra- 
dius of primary; r s relative radius of secondary; k = r s /r p ; J 
central surface brightness ratio; L luminosity; h third light frac- 
tion; e eccentricity; co longitude of periastron. 

3.5.1. V1 8 photometric elements 

Our first photometric solutions for V18 showed that k is not well 
constrained by the light curve without additional external con- 
straints. This is typical for a system with a slightly eccentric 
orbit without a total eclipse. We therefore found solutions for 
the V band light curve which simultaneously fit the light curve 
and spectroscopic values of e sin a> and the light ratio (Sect. I4.ll 
and !4.2l ). including their errors. We could not repeat this proce- 
dure for the R band, since we have no spectroscopic measure- 
ment of the light ratio corresponding to this band. Therefore we 
adopt the V band solution for our final measurements. To get the 
R band surface-brightness ratio and luminosity ratio, we fit the 
R band light curve with the other parameters fixed from the V 
band solution. Table[3]shows our adopted photometric elements. 
Errors are determined using JKTEBOP Monte Carlo simula- 
tions, since residual shift errors indicated that correlated errors 
are not present at a significant level in the obs erved light cu rves. 
Square root limb-darkening coefficients from IClaretl d2000l) for 
our spectroscopically measured metallicity and effective temper- 
atures were used as a starting point. We followed the previously 
mentioned strategy of fixing either the linear or non-linear coeffi- 
cient and solving for the other. The resulting elements remained 
identical to well below 0.01 % whether we solved for the linear or 
non-linear coefficient. Further tests showed that fixing both the 
linear and non-linear limb-darkening coefficients instead of em- 
ploying the spectroscopic e sin a> constraint changes the relative 
radii by only 0.13% and 0.04% for the primary and secondary 
component respectively. This seems to suggest that at least in 
this case the theoretical limb-darkening coefficients are reliable. 
As seen in Table[3] the relative radii of V18 have been measured 
with errors of 0.6% and 0.9% for the primary and secondary 
components respectively. 

3.5.2. V20 photometric elements 

V20 was analysed by GCH08. However, we measured more di- 
rectly and accurately the third light and found it to be outside 
their error estimate. We therefore performed a new analysis us- 
ing their light curve. Table [4] shows our measured photomet- 
ric elements for each band and our adopted solution. From the 
JKTEBOP Monte Carlo and residual shift solutions we found 
that significant correlated errors are present in the observed light 



Table 3. Photometric solution for VI 8 from the JKTEBOP 
code. 



Parameter 


Value 


Spectroscopic constraints 




LJL P (V) 


0.636 ± 0.020 


e sin to 


0.0100 ±0.0011 


Measured parameters 




i O 


89.377 ± 0.022 


r p + r s 


0.0560 ± 0.0002 


k = rjr p 


0.882 + 0.011 


r P 


0.02975 ±0.00018 


r s 


0.02623 ± 0.00024 


e 


0.0193 ± 0.0006 


CO 


149.1 ±2.8 




0.817 ±0.018 


JJJpiR) 


0.791 ±0.074 


LJL„(V) 


0.635 ±0.016 


LJL P (R) 


0.674 ± 0.042 


a (V-mmag.) 


4.36 



curves, and we chose to adopt the larger residual shift errors as 
our error estimates. For the final parameters we adopted an er- 
ror weighted mean of the light curve solutions from each band. 
Due to the systematics present, we do not allow the error of the 
adopted solution to be reduced by averaging the solutions from 
the two bands, but adopt the error from the V band as our final 
estimate. As seen in Table |4] relative radii of V20 have been 
measured with errors of 0.9% and 0.7% for the primary and 
secondary components respectively. Our measurement of r s = 
0.025 19 + 0.00017 is two sigma larger than r s = 0.0248+0.0002 
found by GCH08. Although we also use a different limb darken- 
ing law than GCH08, this has minor effects on the solution. Our 
new measurement of the third light is the main reason we find a 
larger radius for the secondary. 

3.5.3. V80 photometric elements 

As mentioned earlier, V80 shows signs of magnetic activity in 
the light curve and more observations are needed for very ac- 
curate measurements. In any case, to obtain a measurement of 
\ogg to be used for T e $ and [Fe/H] measurements, we did a 
preliminary light curve analysis for V80. To be specific, we first 
rectified the light curves. In Fig. |4]both primary and secondary 
eclipses are covered by observations from two different nights. 
As seen the overall magnitude level remained identical for the 
nights covering the primary eclipse while it changed between the 
two nights covering the secondary eclipse. However, the depth of 
the secondary eclipse did not change between nights. The recti- 
fication we employed was therefore to shift the one secondary 
eclipse coverage which was offset relative to the rest, so that all 
eclipse observations were aligned in magnitude. Then deviating 
out-of-eclipse observations were also shifted to match this mag- 
nitude level. We then found JKTEBOP solutions employing the 
V band light ratio determined from the CMD in Sect. 13.31 The 
results for the relative radii and inclination are shown in Table[5] 
Errors are based on solutions of the V band light curve varying 
the light ratio within the errors given by the external constraint 
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Table 4. Photometric solutions for V20 from the JKTEBOP 
code. 



Parameter 


V 


/ 


adopted 


i (°) 


89.99 


89.99 


89.99 


+0.01 


+0.01 


+0.01 




-0.15 


-0.18 


-0.15 


r P + r s 


0.0700 


0.0710 


0.0702 




±0.0005 


±0.0010 


±0.0005 


k 


0.559 


0.559 


0.560 




i n nnc 
±0.005 


i n nno 
±0.008 


i n nnc 
±0.005 


r P 


0.04491 


0.04551 


0.04501 




i n nnri/in 
±0.00040 


i n nnmo 
±0.000 la 


±0.00040 


r„ 


0.02513 


0.02546 


0.02519 




i n nnn 1 n 
±0.0001 / 


+U.UUUZD 


i n nnn 1 n 
±0.0001 / 


JJJ„(V) 


0.414 


(...) 


0.410 


_i_n run 


(...) 


_i_n ozn 


J si J p\l) 


(...) 


U.J71 


U.J7U 




(...) 


±0.074 


±0.074 


L S /L P (V) 


0.134 


(...) 


0.134 




±0.003 


(...) 


±0.003 


LJL P U) 


(...) 


0.181 


0.181 




(...) 


±0.005 


±0.005 


h{V) 


0.182 


(...) 


0.182 


h{I) 


(...) 


0.197 


0.196 


a (mmag.) 


4.6 


6.8 


(...) 



Table 5. Photometric solution for V80 from the JKTEBOP 
code. 



Parameter 


Value 


Constraints from CMD 




LJL P (V) 


0.15 ±0.05 


Measured parameters 




i (°) 


84 ± 1 


r p 


0.0900 ± 0.0054 




0.061 ±0.012 



and on the level of consistency between V and R band solutions. 
As expected, the errors are much higher than for the other sys- 
tems, in the range of 6-20%. 

4. Spectroscopy 

The spectroscopic observations were carried out in service mode 
with UVES at the ESO VLT during allocation periods 75+77 
(V20) and 81 (VI 8 and V80). Since NGC 6791 is at declination 
+37°, it can be observed at Paranal only at an airmass higher than 
2.1. Therefore all observations were carried out during the short 
observing window of a few hours around meridian passage. Due 
to the faintness of the stars, and in order to minimize slit losses, 
a slit of 1'.'20 width, corresponding to a resolution of approxi- 
mately 37 000, was used. For the observations of V18 and V20, 
the slit was aligned along the parallactic angle (ELEV mode), 
and the atmospheric dispersion corrector (ADC) was not inserted 
in the beam, since it causes a slight loss of flux. The procedure 
for V80 was different, since a fixed slit position allowed us to put 
an additional star on the slit together with V80, to be used for a 
metallicity measurement. The UVES standard 580nm setup, and 



on-chip binning of 2x2 pixels, was used for all observations. The 
corresponding wavelength ranges covered by the two CCD de- 
tectors employed in UVES are approximately 4775-5750A and 
5875-6830A, respectively. The typical S/N per pixel for the red 
chip was between 15 and 25. A total of 13, 15, and 10 usable 
epochs were obtained for VI 8, V20 and V80, respectively, each 
followed by a ThAr exposure for wavelength calibration. 

We used the UVES pipeline to reduce the spectra of all three 
binaries, with a few changes relative to the standard pipeline 
settings: We reduced the threshold to discard deviating pixels 
from 10 sigma to 3 sigma, and we used the semi-automatic flux- 
weighted wavelength calibration. For the observations with more 
than one star on the slit, we applied offsets and reduced the ex- 
traction slit length, to extract the spectrum of each star sepa- 
rately. We did not merge orders, since the regions of overlap 
have low S/N and could introduce unwanted noise in the later 
analysis. 

V20 was previously analysed by GCH08 using the UVES 
pipeline reduced spectra as supplied by ESO. However, inspect- 
ing the raw spectra of V20 revealed that the choice of the ELEV 
mode had the unfortunate effect of allowing light from a close- 
by star (in addition to the third light) to reach the slit for some 
observations. This affects the radial velocity zero-point determi- 
nation significantly, regardless of whether this is done using the 
third light component as in GCH08 or with our present method 
as we describe below. Therefore, these spectra were re-reduced 
with offsets and reduced slit length, to mask out the unwanted 
star. 

The orbital periods of VI 8 and V80 were not well known 
prior to the spectroscopic observations. Therefore, we retrieved 
the data from the ESO Archive the day after each new obser- 
vation, reduced the spectrum, made an orbital solution to find a 
better period estimate, and then rescheduled the remaining ob- 
servations to optimal epochs. 

4. 1 . Radial velocities and spectroscopic elements 

For the radial velocity measur ements we emp l oyed the bro aden- 
ing function (BF) formalism (iRucinskill 19991 l2002i |2004|) . The 
BF method assumes that all components in the observed spec- 
trum can be well described by the same synthetic template spec- 
trum with respect to the relative strengths of spectral lines. This 
assumption holds as long as the spectral types of the compo- 
nents do not differ too much. We used synthetic spectra from 
the grid of lCoelho et alj (120051) and repeated the procedure with 
templates covering a reasonable range in logg, and [Fe/H] 
and found that the effect of the adopted template is much less 
than the final adopted error estimates. We also confirmed the 
similarity of the spectral types of the components by separating 
the spectra (see Sect. 14.31 1. 

In our implementation BFs are calculated for each order and 
averaged to give the final BF profile. This profile is smoothed 
by convolution with a Gaussian. The width of this smoothing 
Gaussian was chosen to minimize the order-to-order scatter of 
radial velocities measured from single-order BFs. Fig. [5] shows 
example BFs for each system. 

The final BFs were then fitted with a Gaussian around each 
peak to measure the radial velocities. We found that this pro- 
cedure gives a lower order-to-order radial velocity scatter than 
fitting a rotational profile, when applied to the individual orders. 
We did not use all spectral orders, since in the bluest region, the 
S/N was too low to obtain a decent BF, and some orders con- 
tain wide lines or telluric lines, and are therefore not well-suited 
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Fig. 5. Averaged broadening functions for the three binaries 
using the optimal templates and Gaussian smoothing, showing 
the characteristics of each system. V80, V20, and V18 BFs are 
shown from top to bottom with offsets of 0.13, 0.06, and 0.00 
for clarity. V18 and V20 are seen to be slowly rotating systems, 
while V80 rotates faster. V20 contains a third light component. 



for radial velocity measurements. We ended up using ten orders 
from each chip. 

Spectroscopic elements were derived from the measured ra- 
dial velocities using the method of L ehman-Filhes implemented 
in the SBOP program (Etzel 2004). The orbital periods were 
fixed at their photometrically determined values. 

Our first solutions had (O - C) errors of up to 1 km/s in 
a systematic pattern identical for the primary and secondary 
components. This indicated that the radial velocity errors were 
dominated by radial velocity zero-point offsets between epochs. 
Since our spectra are calibrated with ThAr spectra taken imme- 
diately after the stellar spectra, these offsets cannot be attributed 
to wavelength calibration errors. Instead they occur because the 
starlight is not exactly centred across the rather wide spectro- 
graph slit, effectively causing an incident angle, which shifts the 
spectrum in wavelength relative to a star perfectly centred across 
the slit. To correct for this shift, we did the following: We iden- 
tified a spectral order with very strong telluric absorption lines 
(A = 6276 - 6299 A). In the adjacent orders we measured the BF 
profile, which is nearly identical for close-by orders. We then 
convolved our synthetic stellar template spectrum, covering the 
wavelength region with telluric absorption lines, with the ob- 
tained BF, to produce the expected stellar spectrum in that order. 
This spectrum was then multiplied with synthetic telluric spec- 
tra which were broadened to the spectrograph resolution, shifted 
in radial velocity from -3 to +3 km/s in steps of 50 m/s, and 
multiplied with factors between 0.5 and 1 .5 in steps of 0. 1 . The 
resulting set of artificial spectra of stellar+telluric lines were all 
cross-correlated with the observed spectra, and the one giving 
the highest cross correlation determined the shift of the telluric 
lines. This shift is our radial velocity zero-point correction, since 
the telluric absorption lines are always at zero radial velocity 
(except for small shifts due to winds in the Earth's atmosphere), 
and the observed shift can only be caused by the aforementioned 
imperfect slit centring. 

Applying the zero-point radial velocity shifts reduced the er- 
rors on the orbital parameters very significantly. It also changed 
the orbital parameters to just outside their original one-sigma er- 
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Fig. 6. Spectroscopic double-lined orbital solution for VI 8. 
Phase 0.0 corresponds to central primary eclipse. The horizontal 
dotted line (upper panel) represents the centre-of-mass velocity. 



ror bars, indicating that these were underestimated because the 
zero-point shifts had not been taken into account. 

In order to check that the errors on the spectroscopic ele- 
ments given by the SBOP program for our zero-point corrected 
radial velocities can be trusted we developed a Monte Carlo rou- 
tine for alternative error calculations. We assumed that the radial 
velocity errors have two Gaussian components, one of which 
represents any remaining radial velocity zero-point error and is 
identical for both components. We subtracted in quadrature the 
standard deviation radial velocity errors found from order-to- 
order scatter from the SBOP (O - C) standard errors in order 
to estimate the magnitude of the radial velocity zero-point error. 
We found this error contribution to be ~ 150 m/s. We then ran 
SBOP a large number of times adding to each radial velocity 
two Gaussian errors with standard deviations of 150 m/s and the 
order-to-order standard deviations, respectively. From the out- 
put solutions we calculated the one-sigma error. We found that 
this approach gives errors consistent with the errors from SBOP, 
and concluded that these errors can be trusted. Adopting a con- 
servative approach we have anyway chosen to adopt the Monte 
Carlo errors, which are slightly (1-10%) larger than the SBOP er- 
rors. F or more details on the error estimates we refer to Brogaard 
J20T0I) . 

We derived orbital parameters for the three binaries from 
both single- and double-lined solutions (meaning that a solu- 
tion was found both by fitting one component at a time, and 
by fitting the two components in a combined solution). For the 
double-lined solutions, the radial velocities of the secondary 
were adjusted to account for the difference in gravitational red- 
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Fig. 7. Spectroscopic double-lined orbital solution for V20. 
Phase 0.0 corresponds to central primary eclipse. The horizontal 
dotted line (upper panel) represents the centre-of-mass velocity. 
Note the different scales in the (O - C) panels. 



Fig. 8. Spectroscopic double-lined orbital solution for the V80. 
Phase 0.0 corresponds to central primary eclipse. The horizontal 
dotted line (upper panel) represents the centre-of-mass velocity. 
Note the different scales in the (O - C) panels. 



shift between components. These adjustments were calculated 
iteratively from the obtained absolute masses and radii and are 
-35 m/s for V18, -180 m/s for V20 and -150 m/s for V80. 
Furthermore, the relative radial velocity weights between com- 
ponents were set according to the difference in (O - C) errors 
between components from single-lined solutions. For V18 the 
orbit is slightly eccentric, so to obtain a consistent set of semi- 
amplitudes corresponding to the same e and at we adopted the 
double-lined solution. For this system, a measurement of e sin u>, 
including its error, was desirable for constraining the light curve 
solution, where this quantity is poorly determined. However, 
since e and oj are correlated parameters, one cannot obtain the 
error on e sin lo from the errors on e and to. We therefore ran 
SBOP using Sterne's method, which fits for e sin w and ecosw 
directly. Sterne's method only works for one component, so 
we did a solution for both the primary and secondary compo- 
nent. This resulted in values for e and a> on either side of the 
adopted solution. We therefore calculated the value of esinw 
from the adopted e and u> while we obtained the error on e sin to 
from the SBOP solutions using Sterne's method. This result is 
e sin u) = 0.0100 ± 0.0011. 

For V20 we cannot rule out a low eccentricity of ~ 0.002 
from the spectra alone, but both V and / light curves prefer an 
orientation of 90 degrees when any eccentricity is forced upon 
them, indicating a strong preference for a circular orbit, which 
is what we adopt. Furthermore, we adopt the single-lined so- 
lutions since they correspond to a consistent set for a circular 
orbit, and this avoids any speculations about differences in ob- 



served system velocity between components due to exact dif- 
ferences in gravitational redshift, convective blueshifts and any 
template mismatch. For V18 such differences must be very small 
because the stars are of very similar masses, radii and spectral 
types, and therefore negligible errors are introduced by adopt- 
ing the double-lined solution for this system. For V80 rotational 
broadening reduces the accuracy of both radial velocities and 
zero-point measurements, but because of the higher radial veloc- 
ity amplitudes caused by the shorter orbital period, the minimum 
masses are still determined to an accuracy better than 1%. 

We list the final spectroscopic elements for the systems in 
Tableland show double-lined solutions in Figs.[6f{8] Tables 15- 
17 containing radial velocities and (O-C) errors can be found on 
CDS. As seen, minimum masses accurate to well below 1 % have 
been obtained for all systems. For V20, the errors on both the 
semi-amplitudes and individual measurements are about half of 
those found by GCH08. Perhaps the most convincing evidence, 
which indicate that our measurements are more accurate, is to 
compare our Fig. Qto their Fig. 8; note in particular the sets of 
measurements close to phases 0.2, 0.25, 0.35, and 0.6, which 
were taken immediately after each other and are therefore ex- 
pected to give identical results. 

As also seen from Table [6] the systemic velocities, y, of the 
systems are in excellent agreement with the value of -47.1 +0.8 
kmjT 1 as determined from 15 cluster members bv lCarraro et al.l 
(2006), who also found the dispersion in the radial velocities of 
the 15 stars to be 2.2 + 0.4 kms~'. This provides a strong argu- 
ment that all systems are cluster members. V18 and V80 are in 
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addition proper-motion members (96% and 86% probability ac- 
cordi ng to Dr. Kyle Cudw orth as stated bv lRucinski et al.ll 1 9961 
and lde Marchi et al]|2007l respectively). Given also the similar- 
ity in the metallicities that we derive (Sect. 14. 3b there should be 
no doubt that these stars are cluster members. 



4.2. Spectroscopic light ratio of V18 

For V18 we found that a spectroscopic light ratio could improve 
the photometric elements (cf. Sect. I3.5.U . For this system we 
therefore found an additional use for our calculated broaden- 
ing functions, since the BF is a proxy for the mean line profiles 
in the stars. The difference in equivalent widths between each 
component peak therefore measures the light ratio. We selected 
BFs from 10 orders for which the mean luminosity ratio corre- 
sponds very closely to the luminosity ratio in the V band accord- 
ing to tests with both Planck functions and synthetic spectra of 
Coelho et al. (2005). For each epoch of observations we fitted 
these BFs with a rotational profile convolved with a Gaussian 
and integrated the area under each peak to measure the ratio of 
equivalent widths between components. If the stars had identi- 
cal spectral types this ratio measures directly the light ratio in 
the V band. Since the light curve solutions suggest an effec- 
tive temperature difference of about 200 K, a very small cor- 
rection of -0.06 was applied to the measured light ratio. We 
found L s /L p = 0.636 + 0.020 where the error is a sum of the 
epoch-to-epoch measurement scatter (+0.005), a contribution 
due to possible differences between light ratio in the V band and 
the spectral region measured (+0.005), an estimated uncertainty 
in the correction due to different spectral types (±0.005), and a 
contribution to account for the exact continuum placement in the 
measurement and the fact that real stars may not be accurately 
described by a rotational profile (+0.005). We chose all these 
error contributions conservatively and add them directly (not in 
quadrature) because they have a systematic character with ex- 
ception of the measurement scatter. We were unable to repeat 
this procedure for the R band, since we found no combination of 
orders for which the light ratio corresponds to the R band light 
ratio. 



4.3. Spectroscopic analysis: T eS and [Fe/H] 

To determine effective temperatures and metallicities of the bi- 
nary stars, we first separ ated their spectra. This wa s done fol- 
lowing the description of iGonzalez & Levatol ((2006), which we 
expanded to handle systems with three components (needed for 
V20). When separating spectra of binary stars the light contri- 
bution from each component to the continuum level cannot be 
determined from the spectra themselves (except when the spec- 
tral types are almost identical, where one can follow the pro- 
cedure we used for V18 in the previous section). We therefore 
employed our photometric light ratios from Tables [3}j5] to cor- 
rect each component spectrum for the continuum contribution 
from the other component. Since these light ratios are measured 
in the Johnson V band, we used Planck functions to calculate 
the corresponding light ratio in each spectral order from initial 
photometric temperature estimates. This procedure was iterated 
with the derived temperatures, with negligible differences. We 
checked that the Planck approximation was adequate, e.g. using 
MARCS models instead produced identical results. The binary 
target V80 was observed with an additional star on the slit on 
some epochs. After acquiring the spectra, a more careful inspec- 
tion revealed that this was in fact two very close stars, which we 



also analyse below. For simplicity we shall refer to this compos- 
ite spectrum as the subgiant (SG) star spectrum. 

We made a detailed spectral analysis of four of the separated 
binary spectra, the primary and secondary of V18 (abbreviated 
V18p and VI 8s), V20p, V80p, and SG. We did not analyse V20s 
and V80s, since the signal-to-noise (S/N) was very low, and 
therefore the placement of the continuum introduced very large 
systematic errors. The continuum was normalized using a syn- 
thetic template spectrum with the assumed atmospheric parame- 
ters of the star. This was done by identifying continuum window s 
with the rainbow software as described in Bruntt et al. (2010b). 

In Fig. [9] we show a small section of the spectra, which shows 
their range of quality. The full-length one-dimensional spectra 
can be found in Tables 18-22 on CDS. The S/N in the continuum 
ranges from 30 (V18s) to 60 (V20p). In Fig. |9]we also overplot 
a high resolution (R=l 10000) high S/N spectrum of a Cen A as 
observed with HARPS. This star has parameters very similar to 
the observed stars in NGC 6791: lBruntt et all d2010al) used the 
HARPS spectrum to determine v sin i = 1.9 + 0.6 km/s, = 
5745 + 80 K, logg = 4.31 ± 0.08, and [Fe/H] = +0.22 + 0.07. 

We used the s emi-automatic soft ware package VWA 
dBruntt et al.l 1201 Oat IClausen et al.1 120081) to analyse the spec- 
tra. The progr am interpolates in a grid of MARCS model at- 
mospheres (Gusta fsson et al.l 20081) and uses atomic line data 
from VALD (Ku pka et al. 19991) . The spectral lines are fit- 
ted iteratively by calcula ting synthetic profiles using SYNTH 
(IValenti & Piskunovl 1 19961) . Oscillator strengths are corrected 
using a s olar spectrum. The software is described more thor- 
oughly in Brun tt et al.l d2010albl) . and we will mention here only 
differences to the standard approach. Most importantly, we made 
the analysis by fixing the log g values as constrained by the ac- 
curate measurements of the eclipsing binaries (Table [101. This 
breaks the degeneracy between log g and T e ff, which is usually a 
problem when analysing spectra at this S/N level, and allows a 
reliable T e s determination. 

As discussed above, the SG spectrum is a composite of the 
spectra of two stars with nearly identical radial velocity. We 
identified them to be stars 1621 and 1627 in the photometry 
of SBG03 (finding chart in Fig. IA.3l l. but their measured mag- 
nitudes are wrong due to their small separation. We therefore 
transformed their magnitudes from our HST photometry into B 
and V magnitudes from which we measured their V light ratio 
and extrapolated along an isochrone to obtain their logg val- 
ues (Table |7J. Since they are very close to the binaries in the 
CMD, this should give logg values with errors of less than 0.05 
dex. The spectrum of SG was then analysed using the a pproach 
for composite spectra described in lClausen et al.l ( 20081) . As for 
the binaries our prior knowledge of log g of the two components 
made it possible to constrain their temperatures. 

We adjusted the vsin; of the synthetic spectra to provide 
the best fit to the observations. The value s of vsin; a r e listed 
in Table [7] We used the calibration from Bruntt et al] (1201 Oal) 
to fix the macroturbulence at 2.0 km/s for all stars. Changing 
the macroturbulence within reasonable errors does not change 
the measured equivalent widths, and thus does not affect the re- 
sults of the spectroscopic analysis. The microturbulence could 
only be determined with rather large errors. However, the values 
found are consistent with the calibration for solar-type stars from 
iBrunttet al.ld2010al) . 

For all the analysed spectra, we adjusted the T e ff of the at- 
mosphere model to make the mean abundances of Fe i and Fe n 
lines agree. We evaluated the internal error on T e ff by changing 
its value until the Fe i and Fe n deviated by one sigma. 
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Table 6. Spectroscopic orbital solutions. 



Parameter 


V18 


V20 


V80 


rixed parameters; 










JtUJ 1 .4JUU 


JJ 1 J 1 -UUU 1 


JtUJZ.JUtJ 


P (days) 


18.798638 


14.469918 


4.88594 


Pree parameters! 








K p (kms -1 ) 


48.094 ± 0.075 


46 900 + 050 


66 05 + 16 




51 523 + 076 


61 <;o + o 10 


87 38 + 26 


7 (km s -1 ) 


-47.743 ± 0.038 


-46.639 ± 0.035 


-45.355 ± 0.098 


e 


0.0196 ± 0.0012 


O.OO(fixed) 


O.OO(fixed) 




149.4 ±3.1 


90.00(fixed) 


90.00(fixed) 


T~)privpH tiaramptprs' 








iWp Sill L ^1V10 ) 


n qck^ -t- n oht^ 


1 nsiss -t- n nrnQ 


i i ^ -t- n nn^Q 


M s sin 3 i (Mq) 


0.9291 ± 0.0032 


0.8276 ± 0.0022 


0.7872 ± 0.0043 


q 


0.9334 ± 0.0020 


0.7615 ±0.0015 


0.7559 ± 0.0028 


a sin (' (Rq) 


37.010 ± 0.040 


3 1.030 ±0.032 


14.819 ±0.029 


Other quantities: 








<Tp(km s" 1 ) 


0.24 


0.14 


0.33 


o",(kms _1 ) 


0.23 


0.31 


0.69 




13 


15 


10 


Time span (days) 


88 


377 


92 



Notes. T is the time of central primary eclipse. 



To confirm the method, iBruntt et"aLl (1201 Oal) compared the 
r e ffS from VWA with results for 10 nearby stars where inter- 
ferometric methods can be used, and found very good agree- 
ment. They claimed a systematic offset of r e ff ; vwA - 7eff;interf. = 
40 + 20 K should be applied to the spectroscopic effective tem- 
peratu res. However comparing the (b - y) temperature calibra- 
tion of Casagr ande et al.l (1201 Ol) for 17 stars in common with the 
spectroscopic analysis of IBruntt et al.1 d2010al) . we find that the 
situation is more complicated. The derived temperature offset is 
sensitive to the absolute magnitude of Vega, which is difficult to 
determine (see Casagran de et alJfeOlOl for a discussion). Equally 
importantly, there are indications that one cannot interpret the 
T e g difference between VWA and other methods as a simple 
offset, since differences may depend on both and [Fe/H]. 
Unfortunately, we do not have enough stars measured to reach 
a firm conclusion. Due to these complications, any attempt to 
apply an offset to the T e ff derived using VWA may introduce 
more systematic error. Therefore, we adopt the T e ff values as de- 
termined from VWA analysis without any correction and add in 
quadrature a systematic error of 70 K to the internal VWA errors. 

In Table [7] we list the adopted log g, the determined T e s , 
vsin/, microturbulent velocity, and [Fe/H] for all stars analysed. 
The errors for [Fe/H] are determined by changing the and 
microturbulent velocity by one sigma in the analysis. The quoted 
errors in this table are internal errors. 

With the final atmospheric parameters determined, we com- 
puted the abundances of all elements with at least four unblended 
lines available. In Table[8]we list the mean abundances and num- 
ber of lines used for each spectrum. We only used relatively 
weak lines with equivalent widths between 10 and 100 mA. The 
quoted errors are the rms of the abundances for each element in 
a given ionization stage. All elements are consistent within 0.1 
dex, except for vanadium, which has a higher mean abundance 
in all stars. 

We tested the effe ct of the assumed hel ium content, follow- 
ing the prescription of lStromgren et al.l (Il982l) . They argue that a 
change in the helium content is effectively the same as a change 
in the pressure, and hence log g, in the atmosphere. Using their 



prescription, we find that even a large change in Y from solar to 
0.33 (high value for NGC 6791) is accommodated by a change 
in logg of only 0.03 which leads to a change in [Fe/H] of 0.01 
or less. At the same time, helium diffusion would tend to lower 
the helium content in stellar atmospheres with time. Thus, for 
the old stars in NGC 6791 more helium will have left the atmo- 
spheres due to diffusion compared to the Sun, and therefore the 
correction due to a different helium abundance should be even 
smaller and effectively negligible. Finally, we repeated the anal- 
ysis of the stars when changing the light ratio by two sigma, and 
found negligible changes to the determined T e ff and metallicity. 

For our final estimate of [Fe /H] for NGC 679 1 we calculated 
the weighted mean of the six stars [Fe/H] in Table[7] resulting in 
[Fe/H] NG c 6791 = +0.29 + 0.03 + 0.07, where the quoted uncer- 
tainties are the weighted mean error and the adopted systematic 
error. 

There are several spectroscopic measurements of the metal- 
licity of NGC 6791 in the literature, with results spanning a 
rather wide range in [Fe /H] from +0.30 to +0.47. Nearly all are 
all ba s ed on giant stars jFriel & Janeslll993l iPeterson & Greeril 
1998; iGratton et all 120061: lOrigha et al.l l2006t ICarraro etaLl 
2006). These analyses of giant stars are complicated by the 
presence of molecules and blending of lines, especially since 
they are based on low-to-medium resolution spectra (all have 
R < 30 000). In such cases, blen ding makes the plac ement of the 
continuum difficult ( cf. Fig. 2 inlGratton et alj200 6 and Fig. 2 in 
ICarraro et al.l l2006). Boesgaard et al. (2009) analysed two sub- 
giant stars at higher resolution (R = 45 000), but they adopted a 
log g which is 0.3 dex higher than inferred from our binary mea- 
surements. Interestingly, their different log g and T e g parameter s 
conspire to give [Fe/H] = +0.30 ± 0.08 dBoesgaard et al.ll2 009). 
in agreement with our result, despite their different atmospheric 
parameters. 

We believe that the metallicity determined from our sample 
of stars is reliable, since the parameters of the stars are close to 
the Sun, and therefore less prone to systematic effects coming 
from the adopted model atmospheres. Furthermore, our use of 
accurate log g values from the binary measurements has strongly 
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constrained one free parameter in the spectroscopic analysis. We 
note also that the [Fe/H] measurements of the two single stars 
in the SG spectrum give results in excellent agreement with the 
binary stars. This is clear evidence that the spectral separation of 
the binary stars has not affected the metallicity measurements. 

5. Absolute dimensions and distance 

5.1. The reddening of NGC 6791 

The reddening, E(B - V), is needed to obtain the true distance of 
the cluster, and to constrain isochrones in the CMD. Using V18p, 
VI 8s, and V20p, we restricted the acceptable range of E(B - V) 
by requiring our spectroscopic values of T e ff to agree with the 
temperatures obtained from a photometric c olour-temperature 
calibration. We used the calibrations from ICasagrande et al.l 
(201( j) and from MARCS model atmospheres (calculated for use 
inlVandenBerg et al.l2010l) . For a range of E(B-V) we found the 
de-reddened colours of the binary components, and from those 
the photometric temperatures, assuming [Fe/H] = +0.3. The two 
calibrations give very consistent T e ffS with the largest deviation 
being 40 K. We set the one-sigma limits on E(B - V) where 
differences between the spectroscopic and the photometric r e ffS 
become greater than 200 K for one of the stars, or greater than 
100 K for all three stars. This one sigma range in E(B - V) is 
0.135-0.185 with the best match at 0.16. Repeating the proce- 
dure for an assumed [Fe/H] of +0.2 and +0.4 showed that the 
result is only slightly dependent on [Fe/H] and thus the accept- 
able range in E(B - V) is independent of [Fe/H] within +0.1 
dex of our [Fe/H] measurement. We thus find from this method 
that E(B — V) = 0.160 ± 0.025 in excellent agreement with other 
recent results (VandenBerg et al. 2010; and references therein), 
and our other measurement, which we now describe. For an ad- 
ditional measure of the inter stellar reddening, we used the cali- 
bration of iMunari & Zwitter] (1 19971) for the relationship between 
the equivalent width (W) of interstellar Nal D lines and E(B-V). 
We first combined all our binary spectra at their barycentric ve- 
locity in order to average out the spectral signal from the stars. 
In the combined spectrum the continuum has been suppressed 
to the left of the interstellar lines because of the Nal D lines 
from the binary stars, which are at different wavelengths in each 
spectrum, because the stars orbit each other. We assumed that 
the continuum suppression above a line across the interstellar 
Nal D line removes the stellar contribution to the interstellar line 
depth. Removing first this contribution, we measured the equiv- 
alent widths of the interstellar Nal D lines. We repeated the pro- 
cedure for different selections of spectra, first all binary spectra, 
then selecting only spectra where deep lines are as far from the 
interstellar line as possible, and then this same procedure using 
only spectra from one star. 

Munari & Zwitter (1997) find that W(D1)/W(D2)> 1.1, 
which enables a minimum estimate of W(D1) from W(D2). 
Since in our case this minimum estimate is slightly higher than 
the direct measurement, we use the mean of the two measure- 
ments as our final measurement. Results are shown in Table [9] 
from which we adopt 0.354+0.010 as the measured W of the Nal 
Dl line giving most weight to the measurement using many bi- 
nary spectra while avoiding strong stellar contamination. Using 
this value with Table 2 of Munari & Zwitter (1997) we find 
E{B - V) = 0.153 ± 0.010. Munari & Zwitter state that for 
E(B - V) below 0.5 the reddening uncertainty is ~ 0.05, but 
from their Fig. 5 it is clear that for E(B - V) below 0.2, the 
uncertainty is even smaller. Despite this we conservatively adopt 
E(B - V) — 0. 153 ± 0.050 as our measurement from this method, 



which is admittedly not very precise but in excellent agreement 
with the above result. 

We adopt the most accurate of our two measurements to be 
used in the following: E(B - V) = 0.160 ± 0.025. 

5.2. Absolute dimensions and distance of V18, V20 and V80 

Absolute dimensions for VI 8, V20, and V80 are calculated from 
the elements given in Tables [3]-[6] and listed in Table |T0j As 
seen, masses and radii have been obtained to an accuracy well 
below 1% for both V18 and V20. Errors on the masses are in 
the range 0.27-0.36%, while errors on the radii are in the range 
0.61-0.92%. These measurements are therefore of an accuracy 
comparable to, or b etter than, that of field star eclipsing binaries 
dTorres et alJl2010l) . The masses of V80 are also measured with 
errors below 1%, but as expected the radii are not very precise. 
However, the error in log g for the primary component is still 
only +0.05, which was useful for the [Fe/H] measurement. Our 
measurements of the radius of the V20 secondary and the mass 
of the V20 primary are outside the 1 cr error estimates of GCH08. 
These differences are due to our improved measurements of the 
third light and radial velocity zero-points. 

As also seen in Table [TOl our measurements of the rotational 
velocities are in agreement with the theoretical synchronous ro- 
tational velocities. V18 has an eccentric orbit, but the difference 
between theoretical synchronous and pseudo-synchronous ve- 
locities is only 0. 1 km/s and has been included in the theoretical 
values by increasing the error of the theoretical prediction. 

From our measured radii and spectroscopic r e ff values, we 
calculated the luminosities of V18p, VI 8s and V20p, which to- 
gether with a bolometric correction (B.C.) gave My of the stars. 
From these, and our measured V and E(B - V), we calculated 
the distance moduli and true distance, which we also tabulate 
in Table \10\ As seen the distance moduli of the three compo- 
nents agree very well, with the error in the T e ff values being the 
dominant error contribution. We used B.C. from MARCS stel- 
lar atmospheres, which are in agreement with empirical B.C.s 
from Flower (1996) within 0.02. The true distance moduli were 
found using our measurement of E(B — V) = 0.160 + 0.025 
and A v = 3.09 x E(B - V) <lMcCaUll2004h . From the indi- 
vidual distance moduli of V18p, VI 8s, and V20p, we calcu- 
lated the apparent and true distance moduli of NGC 6791 to be 
(to - M) v = 13.51 + 0.06 and (V - My) = 13.01 + 0.08. The 
distance is then d(NGC 6791) = 4.01 + 0.14 kpc. 



6. A first comparison with theoretical models 

GCH08 used their measurements of V20 to show that the dif- 
ference in its predicted age due to the specific stellar models 
that were employed was about four times greater than their mea- 
surement precision. However, they were unable to determine 
which of the models (if any) were the most trustworthy ones. The 
present study has significantly improved the measured properties 
of the binaries — including [Fe/H], T^, and E(B - V); conse- 
quently, even tighter constraints may be placed on the models. 
Of particular importance is the fact that, given our accurate mass 
and radius measurements for both V18 and V20, we can define 
four points (and their uncertainties) on the mass-radius diagram 
and therefore constrain the curvature of the isochrone that best 
reproduces the observations. 

Although we could repeat the analysis carried out by 
GCH08, an examination of the models that they used re- 
vealed that the different assumed helium mass-fraction, Y, of 
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Table 7. Parameters determined from the spectral analysis (log g values were fixed). The uncertainties are internal errors. 





V18p 


V18s 




V20p 


V80p 


SGp 




SGs 


T eB [K] 


5600 ± 50 


5430 i 


:95 


5645 ± 50 


5600 ± 95 


5540 i 


:60 


5570 ± 75 


l°gg 


4.35 


4.43 




4.18 


4.21 


4.04 




4.19 


v sin i 


4.00 ± 1.00 


3.50 ± 


1.00 


4.50 ± 1.00 


14.00 ± 2.00 


4.00 ± 


1.00 


4.00 ± 1.00 


^micro 


1.00 ±0.12 


0.95 ± 


0.10 


0.90 ±0.10 


1.10 ±0.10 


1.20 + 


0.12 


1.25 ±0.14 


[Fe/H] 


+0.31+0.06 


+0.22 


±0.10 


+0.26 ± 0.06 


+0.34 ±0.10 


+0.32 


±0.07 


+0.30 ± 0.08 



Table 8. Abundances and number of lines used in the spectral analysis. The uncertainties are the standard deviation of the mean. 





V18p 


V18s 


V20p 


V80p 


SGp 


SGs 


Si i 


+0.35 ± 0.05 


9 


+0.46 


±0.06 


4 


+0.23 ± 0.03 


10 


+0.32 + 0.04 


6 










Til 


+0.26 ± 0.04 


11 


+0.28 


±0.05 


6 


+0.27 + 0.06 


7 






+0.17 ±0.04 


8 


+0.11 ±0.07 


8 


Vi 


+0.48 ± 0.07 


7 


+0.58 


±0.05 


8 


+0.35 ± 0.03 


8 


+0.57 ± 0.07 


5 










Cri 


+0.11 ±0.05 


4 
























Fe i 


+0.31 ± 0.01 


63 


+0.22 


±0.03 


50 


+0.26 + 0.01 


76 


+0.34 + 0.03 


30 


+0.32 ± 0.02 


56 


+0.30 + 0.02 


56 


Fe ii 


+0.31 + 0.04 


5 


+0.21 


±0.04 


5 


+0.26 + 0.08 


6 


+0.34 + 0.05 


5 


+0.31 ±0.02 


6 


+0.31 ±0.03 


6 


Nil 


+0.38 + 0.03 


17 


+0.35 


±0.05 


18 


+0.31 ±0.02 


18 


+0.42 + 0.06 


10 


+0.42 ± 0.05 


10 


+0.34 + 0.04 


10 



6150 



6155 



6160 



6165 




6150 6155 6160 6165 

Wavelength [A] 

Fig. 9. The five spectra analysed with VWA are compared. The continuum level of the spectra is offset in steps of 0.25 for clarity. 
Each spectrum is compared with the HARPS spectrum of a Cen A (thin red line). 



NGC 6791, is the pri mary reason fo r the differences in the 
ages that they derived dBroga ard 2010). Since most of the cur- 
rently available isochrones do not contain variations in Y at fixed 
[Fe/H] values, a meaningful comparison of the models produced 
by different workers cannot be made. Moreover, there have been 
a number of recent revisions to nuclear reaction rates and other 
physics ingredients (see the next section), which further limits 
the usefulness of comparisons with most published models be- 
cause they are not based on the most up-to-date stellar physics. 
In this paper, we have therefore chosen to make some initial 
comparisons with isochrones for a range in [Fe/H] and Y, based 
on the tracks computed using just one stellar evolutionary code 
that has incorporated recent advances in nuclear reactions, opac- 
ity, etc. 

6.1. Stellar models 

The stellar models that are used in the present and forthcom- 
ing analysis of the binary stars V18 and V20 in NGC 6791, to- 
gether with the cluster CMD, were generated using a signifi- 



cantly updated version of t he University of Victoria evolution- 
ary c ode (last described by VandenBerg, Bergb usch. & Dowlerl 
2006; and references provided therein). The most important im- 
provement that has since been made to it is the incorporation 
of the treatment of the gravita tional settling of helium presented 
bv iProffitt & Michaudl (Il99~il) . A simple prescription for turbu- 
lent mixing has also been implemented, given that additional 
mixing below envelope convection zones must be assumed in 
order for theoretical models to reproduce the observed varia- 
tions of the surface metal abundances of stars between the main 
sequence and the giant br anch in old stellar systems (e.g., see 
the study of NGC 6397 bv lKorn et ani2007l) . The free parame- 
ter in the adopted formulation of turbulence is set by monitoring 
the settling and nucleosynthesis of lithium and requiring that a 
Standard Solar Model yield the observed Li abundance of the 
Sun. In fact, it turns out that, with such a normalization, models 
for very metal-poor stars predict variations of the surface he- 
lium abundance with evolutionary state that are quite similar to 
those reported by Korn et al. For a much more detailed descrip- 
tion of these results and the adopted treatment of diffusion and 
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Table 9. Equivalent width measurements of the interstellar Nal D lines. 



Spectra used 


W(D1) 


W(D2) 


W(Dl)min. from W(D2) 


W(D1) from W(D1)+W(D2) 


V18+V80+V20 


0.3406 


0.3200 


0.3520 


0.346 


(V 1 8+ V80+ V20)selected 


0.3544 


0.3222 


0.3544 


0.354 


V80selected 


0.3309 


0.3558 


0.3913 


0.361 



Notes. W(D1) and W(D2) are direct measurements of the Nal lines. W(D1) min. from W(D2) is a minimum value of W(D1) calculated from 
W(D1)/W(D2) > 1.1. W(D1) from W(D1)+W(D2) is the mean value of the two estimates of W(D1). 

Table 10. Astrophysical data for VI 8, V20 and V80. 





VI 


8 


V20 


V80 




Primary 


Secondary 


Primary 


Secondary 


Primary 


Secondary 


Absolute dimensions: 














M/Mq 


0.9955 ± 0.0033 


0.9293 ± 0.0032 


1.0868 ± 0.0039 


0.8276 ± 0.0022 


1.0588 ±0.0091 


0.8003 ± 0.0062 


R/Rq 


1.1011 +0.0068 


0.9708 ± 0.0089 


1.397 ±0.013 


0.7813 ±0.0053 


1.341 ±0.081 


0.90 ±0.18 


l°gg (cgs) 


4.3524 ± 0.0053 


4.4319 ±0.0080 


4.1840 ±0.0078 


4.5698 ± 0.0059 


4.208 ± 0.052 


4.43 ±0.18 


T eB (K) fl 


5600 ± 95 


5430 ± 125 


5645 ± 95 


4824 ± 140 


5600 ± 95 


(...) 


v rot (kms -1 )* 


4.0+1.0 


3.5 ± 1.0 


4.50 ± 1.00 


(...) 


14.1 ±2.0 


(...) 


Vsync (kms-'f 


2.96 ±0.14 


2.61 ±0.13 


4.88 ± 0.04 


2.73 ± 0.02 


13.88 ±0.83 


9.33 ± 1.9 


V 


18.247 ±0.020 


18.740 ± 0.020 


17.695 ±0.020 


19.874 ±0.020 


(-Y 


(...y 


M v 


4.76 ±0.10 


5.22 ±0.10 


4.19 ±0.10 


6.44 ±0.18 


(-Y 


(...y 


(V-M v ) 


13.49 ±0.10 


13.52 ±0.10 


13.50 ±0.10 


13.44 ±0.18 




U) c 


(Vo - M v ) 


13.00 ±0.13 


13.03 ±0.13 


13.01 ±0.13 


12.95 ± 0.20 


(...y 


(...y 


distance (kpc) 


3.98 ± 0.25 


4.04 ± 0.25 


4.00 ± 0.25 


3.89 ±0.34 







Notes. 

T eB values are our spectroscopic measurements, except for T eff of V20s, which is calculated as T cff of V20p minus the difference in photometri- 
cally predicted temperatures between V20p and V20s assuming E(B - V)=0. 16. {b) v ro , is the observed equatorial rotational velocities calculated 
by combining vsini and sin; measurements. v sync is the theoretical equatorial velocity for synchronous rotation. (c ' Absolute magnitudes and 
distance moduli have not been calculated for V80 because these are very uncertain due to uncertain V magnitudes and radii (and T c g for V80s) 
caused by the magnetic activity. 



turbulence, we refer to a forthcoming study which examines the 
impact of varying the abundances of individual heavy elements 
(D. VandenBerg et al., in preparation). The settling of heavy ele- 
ments is not considered. However, as discussed below, this omis- 
sion is expected to have only minor consequences for the derived 
age of NGC 6791. 

Since 2006, there have been a number of improvements to 
the rates of several reactions in th e pp-chain a nd the CNO-cycle 
(for a convenient summary, see [Weiss 2008]) — including, in 
particular, the rate of the "bottle-neck" u N(p, y) 15 reaction 
(Martaeta 1. 2008). The current version of the Victoria code 
has been revised to take these advances into acc ount, as well 
as the improved conductive opacities published bv lCassisi et alJ 
(2007). In addition, particular care has been taken concerning 
the surface boundary conditions that are needed to s olve the stel- 
lar structure equations (see, e.g..lVandenber gfl992h . because the 
predicted scale depends so sen sitively on them. As demon- 
strated by 1 VandenBerg et akl d2008l) . simple model atmospheres 
that are obtained by integrating the hydrostatic equation i n con- 
junction with the scaled empirical iHolwege r & Muellei] (1 19741) 
T-t relation yield pressures at T = T e g that are in very good 
agreement with the values obtained from scaled, differentially 
corrected MARCS model atmospheres over wide ranges in T e g , 
log g, and metallicity. That is, the former approach provides a 
good approximation to the use of modern, blanketed model at- 



mospheres as boundary conditions; hence, we have opted to de- 
termine the boundary pressures in t his way, using the fit to th e 
Holweger-Miiller T-t data given bv | Vandenberg & Polll (fl 989). 
In sup port of this choice, VandenBerg, Casagrande, & Stetson 
d2010t) have shown that the resultant temperatures predicted by 
stellar models for field subdwarfs are in excellent agreement 
with those determined for them using the infra-red flux method. 

Fig. [10] shows that the Victoria code produces evolutionary 
tracks for super-metal-rich (SMR) stars that agree quite well 
with those computed for the same masses and initia l helium con- 
tent u sing the completely independent MESA code dPaxton et al.1 
2010; see also http://mesa.sourceforge.net). Except for the 
fact that the MESA code takes the gravitational settling of 
the metals into account (i.e. , as well as helium) using the 
iThoul. Bahcall. & Loebl d 19941) formalism, the physics incorpo- 
rated in the two codes is quite similar. The MESA computations 
assume an initial mass-fraction abundance of the metals, Z, that 
is about 20% higher than the predicted surface metallicity at an 
age of ~ 7 Gyr, when the corresponding [Fe/H] value is w +0.30. 
The surface value of Z in the Victoria models does not vary with 
time. Encouragingly, the tracks superimpose one another nearly 
exactly between the zero-age main sequence and the lower red- 
giant branch (RGB), and they predict the same ages at the tip 
of the RGB to within 1 .4%, despite the noted differences in the 
treatment of diffusion. Thus, models computed using either code 
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Fig. 10. Comparison between Victoria and MESA stellar evo- 
lutionary tracks for masses and chemical abundances (as indi- 
cated) that are relevant to an investigation of NGC 6791 and the 
binary stars V18 and V20. 



Fig. 11. The measured masses and radii of the components of 
V18 and V20 compared to best matching isochrones for [Fe/H] 
= +0.20 and helium contents of Y = 0.27, 0.30, and 0.33. From 
left to right the stars are V20s, VI 8s, V18p, and V20p. Error 
boxes are two-sigma. 



can be expected to predict very similar temperatures and ages 
for the NGC 6791 binaries considered in this study, which are 
all in their core H-burning phases. Although the settling of the 
metals reduces the minimum mass that re tains a convective cor e 
throughout its main-sequence phase (see iMichaud et all 12004). 
NGC 6791 appears to be old enough (see below) that this is 
unlikely to be a concern for th e interpretation of the observa - 
tions. Indeed, as pointed out bv IChristensen-Dalsgaardl (12009). 
the ~ 50% reduction in the rate of the 14 N(p, y) ls O reaction 
(Marta et al. 2008) has the important consequence of increas- 
ing the mass of the lowest mass star that retains a convective 
core throughout its main-sequence phase, thereby decreasing the 
age of the oldest isochrone that predicts a gap near the turnoff. 
The main reason for the difference in the predicted location of 
the RGB is that the MESA code requires a lower value of the 
mixing-length parameter (q-mlt ~ 1 -8) than the Victoria code 
(* 2.0) to satisfy the solar constraint. 

For this initial analysis of NGC 6791 and its binaries, we 
have adopted the relative ab undances of the heavy elements re- 
cently derived for the Sun by Aspl und et alj {2009). 

6.2. Mass-Radius diagrams 

The best starting point for an analysis of binary stars is ar- 
guably the mass-radius (MR) diagram because it provides the 
most direct comparison between theory and observations that 
can be made, being independent of uncertainties in distance, red- 
dening, and colour-temperature transformations. Since the mea- 



sured [Fe/H] value is very close to +0.3 (with a systematic er- 
ror that is difficult to determine), we will compare our obser- 
vations with models that have [Fe/H] = +0.2, +0.3, and +0.4, 
with Y = 0.27, 0.30, and 0.33 at each [Fe/H]. In Figs.[TLHl3]the 
measured masses and radii of the components of V18 and V20 
are compared with isochrones for these chemical abundances (as 
noted) and for the indicated ages, which span the range needed 
to match the most massive, and thus most highly evolved, star 
in the sample (which is the primary of V20). It is worth reiter- 
ating the fact that having four measurements in the MR diagram 
allows us to constrain the shape of the isochrone on this plane. 

The error boxes correspond to two sigma uncertainties. We 
have chosen to do this because, with an increasing number of 
measurements in the mass-radius diagram, the likelihood of fit- 
ting all of the observations to within one sigma decreases. On the 
other hand, simple statistics cannot be used to predict the prob- 
ability of fitting all of the observations to within 1 cr (or 2 <x) 
given that errors in the masses and radii can compensate for one 
another to a extent that depends on the slope of the isochrone 
in the vicinity of each of the observed points. For more rigorous 
results, we intend to use Monte Carlo simulations to quantify 
how well we should expect to fit our observations (Brogaard et 
al. 2010, in prep.; hereafter Paper II). 

Figs. [TTf[T3l show quite clearly that the morphology of an 
isochrone over the range in mass encompassed by the binaries 
depends quite sensitively on the assumed helium abundance (and 
only weakly on [Fe/H]). Indeed, it is only because we have mea- 
surements of multiple eclipsing binaries that we are able to use 
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Fig. 14. Measurements of the detached eclipsing binaries V18 and V20 and the cluster (V - I)q, V CMD compared to isochrones 
with [Fe/ H] - +0.30 and Y = 0.33. The observations in the CMD are from a 2010 re-reduction of the SBG03 data set by P. Stetson 
(seefVandenB erg et alJl20ld Only the stars with the highest precision VI photometry have been plotted. 



the shape of an isochrone on the MR plane to constrain the value 
of Y. This would obviously not have been possible if only VI 8 
(the two middle measurements) or only V20 (the highest and 
lowest points) were used in the comparisons with theory — as 
a change in Y would then be indistinguishable from a change 
in age. With four stars, this degeneracy is now broken. As seen 
from the three figures, isochrones for Y = 0.27 are unacceptable 
for all choices of [Fe/H], because the isochrones miss one star 
by more than two sigma, and only touches the others at or near 
the corners of their 2 <x error boxes. Isochrones for Y = 0.30 
match the observations somewhat better, though they also tend 
to pass through the corners of the 2 cr error boxes: the least ac- 
ceptable of the Y = 0.30 isochrones is that for [Fe/H] = +0.40. 
In fact, the shallow curvature suggested by the measured masses 
and radii indicates a preference for the highest helium abundance 
that we have considered (Y = 0.33), independently of the as- 
sumed [Fe/H] value. However, such a high value of Y requires 
rather young ages, which turn out to be problematic, as we will 
now demonstrate. 

6.3. T ef f-luminosity and colour-magnitude diagrams 

The isochrones plotted in the various MR diagrams that we have 
considered may, of course, also be compared with observations 
of the components of V18 and V20 on the log(r e ff )-log(L/L©) 



(TL) diagram, using luminosities for the latter which can be 
readily calculated from their measured radii and the values of 
T e ff that were derived in our spectroscopic analysis. In addition, 
the same isochrones may be fitted to the CMD of NGC 6791 to 
determine the apparent distance modulus that corresponds to a 
given age. The value of (m - M)y derived in this way should 
agree with the modulus found from the binary stars themselves 
(13.51+0.06, see TablefTOt. An acceptable model must therefore 
be able to fit the MR and TL diagrams, as well as the observed 
CMD, simultaneously. 

Fig.fl4lprovides an example, for the case of [Fe/H] = +0.30 
and Y = 0.33, that can be ruled out because it is not possible 
for the models to match the observations on these three planes 
in a fully consistent way. Clearly both the predicted effective 
temperatures and the apparent distance modulus are way too 
high. For the same reasons, the other Y = 0.33 cases (as well 
as the possibility of even higher helium contents) can be ex- 
cluded. Because our examination of the MR diagram has already 
precluded Y < 0.27 as viable possibilities, we are then left to 
consider the models having Y - 0.30. If [Fe/H] = +0.20, an 
isochrone for the age that does the best job of reproducing the 
MR diagram will fit the turnoff of the cluster CMD only if the 
apparent distance modulus is ~ 13.95 (not shown), which is un- 
acceptably high. As already noted, models for [Fe/H] = +0.40 
provide an unsatisfactory fit to the mass-radius data, leaving us 
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Fig. 15. As Fig.[14]but for isochrones with [Fe/H] = +0.30 and 

with only the Y = 0.30, [Fe/H] = +0.30 case. As illustrated 
in Fig. Q3] these parameter choices result in a level of agreement 
between theory and observations that can be considered satisfac- 
tory. The failure of the isochrones to reproduce the lower main- 
sequence of NGC 6791 on the CMD could well be due mostly 
to problems with the adopted colour-r e ff relations. The discrep- 
ancies would be significantly smaller had we used the empirical 
colour transformations given by Casagrande et al. (2010); see 
VandenBerg et al. (2010). For stars brighter than My ~ 6 the 
MARCS and the Casagrande et al. transformations are in good 
agreement. 

While we could explore the impact of small changes in Y and 
[Fe/H] close to this solution, we should first expand our com- 
parisons with theoretical models to those which assume differ- 
ent chemical abundance patterns. Fo r instance, we have noticed 
that models which assume the older iGrevesse & Sauvall {l998) 
metals mixt ure, which are also fa voured from helioseismological 
studies (e.g jBahcall et al.l2005l) . are able to reproduce the prop- 
erties of the secondary of V20 much better than those that as- 
sume the Asplund et al. (2009) mix of heavy elements. (Whether 
better agreement can also be found for the other binary compo- 
nents remains to be seen.) In addition, it is worthwhile to ex- 
amine the effects of varying the CNO abundances, given that 
they have been found to be und er- abundant in NGC 6791 by 0.2- 
0.3 de x relat ive to scaled solar (ICarretta et al.l2007tlOriglia et al.l 
120061) while lAnthonv-Twarog et al.l (120071) also find indications 
of C being under-abundant. It is well known that the CNO el- 
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ements have important consequences for stars in their core H- 
burning phases through both nucleosynthesis and opacity effects. 
In Paper II we will present a thorough investigation of these and 
other issues, in the hope that we will obtain much more satis- 
factory fits of stellar models to both the CMD of NGC 6791 and 
the cluster binaries than those reported here. Thus, our best es- 
timates of the age and helium abundance of NGC 6791 will be 
reported in Paper II. 

7. Summary and conclusions 

In this paper we have presented extensive photometric and spec- 
troscopic observations of the eclipsing binaries VI 8, V20 and 
V80, which are members of the old open cluster NGC 6791, and 
determined accurate masses and radii for the components of two 
of these systems. Additionally we exploited the eclipsing bina- 
ries for reliable measurements of spectroscopic T e gs, and the 
metallicity and reddening of NGC 679 1 . 

By performing a combined stellar model comparison with 
multiple eclipsing binaries in MR and TL diagrams and the clus- 
ter CMD we showed that we can constrain stellar models better 
than ever before. This allowed us to constrain the helium con- 
tent, and thereby obtain a more precise (and hopefully accurate) 
age of the cluster, although in a model-dependent way. It turned 
out that in order for models to match both the MR and TL di- 
agram and the CMD, only an age close to 7 Gyr is allowed by 
the observations, even though models are selected with a range 



16 



K. Brogaard et al.: Eclipsing binaries in the open cluster NGC 6791 




Ji /Ji 



Fig. 12. As Fig.EJbut for isochrones with [Fe/H] = +0.30 
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Fig. 13. As Fig.[TT]but for isochrones with [Fe/H] = +0.40 



of ±0.1 dex in [Fe/H] and a large range in Y. Our best current 
estimate of the cluster age and helium content is 7.0 Gyr and 
Y=0.30. This solution indicates that the helium enrichment law 
is AY/AZ ~ 2, assuming that such a relation exist and applies to 
NGC 6791. However, these conclusions may depend on the de- 
tails of the stellar models and the detailed abundance pattern of 
NGC 6791, which will be investigated in the forthcoming paper 
II. 

We would like to point out that NGC 679 1 c ontains a large 
number of additional detached eclipsing binarie s (Rucinski et alJ 
119961: Ide Marchi et alj 120071: iMochejska et al.ll2005l ). Accurate 
measurements of these, which we have shown to be possible, 
would allow even tighter model constraints in the MR diagram 
by extending the mass range over which the observed mass- 
radius relation must be reproduced. This has the potential to 
strengthen tests of other physical aspects of stellar models. This 
is so because incorrect model physics may be concealed in the 
models due to the flexibility offered by the uncertain helium con- 
tent, but only up to a certain level. Furthermore, N GC 679 1 is in 
the F ield-Of-View of the NASA Kepler mission dBorucki et al.l 
2010), which will not only allow many more detached eclips- 
ing binaries to be found, and their periods determined, but also 
complementary model constraints from asteroseismology of the 
giant stars in the cluster. 

Multiple detached eclipsing binaries have also been con- 
firmed in a number of other old open clusters (Brogaard 2010; 
Talamantes et al. 2010). Extending this kind of analysis to these 
clusters will allow more reliable cluster ages to be determined 
and more aspects of stellar models to be tested, when stellar 
models have to reproduce observations in clusters with a wide 
range in age and metallicity. By extending studies also to the 
globular clusters, detached eclipsing binaries could ultimately 
provide the strongest constraints on their ages and helium con- 
tents. 
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Appendix A: Finding charts 
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Fig. A.l. Finding chart for VI 8. Plate scale is 0.19"per pixel. Image from ALFOSC at the Nordic Optical Telescope. 
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Fig. A.2. Finding chart for V20. Plate scale is 0.05 "per pixel. Image from ACS on HST 
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Fig. A.3. Finding chart for V80 and SG. Plate scale is 0.05 "per pixel. Image from ACS on HST 
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